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The vacuum gluonic correlation length is extracted from recent lattice data on spin-dependent interquark 
potentials in heavy quarkonia. It is shown that the data are consistent with extremely small values of the 
correlation length, Tg < 0.1 fm. 
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1. INTRODUCTION 

The history of the nonperturbative vacuum gluonic 
fields in QCD can be started with the introduction of 
condensates in the framework of the QCD sum rules 
[1], where the notion of the nonperturbative vacuum 
fields was introduced and the first estimate of the glu- 
onic condensate, Ga = ^{F^^F^^) = 0.012 GeV^, was 
given. Higher condensates and other possible vacuum 
averages of local operators, in line with the Wilson ex- 
pansion, were introduced and estimated in Ref. |2j (for 
a review see Ref. [3]). However, for understanding the 
,nonperturative dynamics of the QCD vacuum another 
characteristic is vitally important — namely, the vac- 
lUum correlation length of gluonic fields Tg, which de- 
fines the nonlocality of gluonic excitations. At a phe- 
nomenological level it was discussed in Ref. [4j, while 
its rigorous definition was given in the framework of the 
Field Correlator Method (FCM)i) [3 [B] (see Ref. [7] 
|for a review of the method). The physical role of Tg 
for the phenomenology of hadrons is quite important: 
,in particular, for hadrons of the spatial size R and the 
temporal size Tq, the QCD sum rule method can be ap- 
plied if R, Tq <SiTg, while potential- type approaches are 
valid in the opposite limit, R,Tq ^ Tg. An example 
of the exactly solvable theory with Tg = is provided 
by the 't Hooft model for QCD in two dimensions [8] 
(see Ref. [S] for a review), which reveals a lot of inter- 
esting phenomena reminiscent of those one expects in 
four-dimensional QCD. 

Direct lattice measurements of the gluonic corre- 
lation length give rather small values for the latter, 
Tg K, 0.15 -^ 0.3 fm [1^. In this paper we adopt another 
strategy: using the FCM we extract the vacuum corre- 
lation length from the spin-dependent potentials mea- 
sured recently on the lattice, in the quenched approxi- 
mation and without cooling [11] . We argue that the data 
prefer even smaller values of the vacuum correlation 
length, Tg < 0.1 fm. Notice that the spin-dependent 



potentials possess a direct physical meaning and are ex- 
pected to be free of any artifacts of the definitions and 
methods used. In particular, the lattice field strength 
correlators, also measured in Ref. [TT], are calculated 
as matrix elements of operators defined in (potential) 
nonrelativistic QCD and simulated on the lattice. They 
are to be multiplied by the appropriate renormalisation 
factor in order to give the corresponding correlators in 
continuum. Although such renormalisation factors are 
expected to be calculated directly from QCD, they are, 
generally speaking, not known at present. 



1. SPIN-DEPENDENT POTENTIALS IN THE 
FIELD CORRELATOR METHOD 

Spin-dependent interaction in a heavy quarkonium 



is well known 
reads: 
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where 777^ and cri (* — 1,2) are the quark masses and 
spins, respectively. Primes denote the derivatives with 
respect to the interquark separation r. The static in- 
terquark potential Vb(r), together with the potentials 
Vi{r) and ^2(7'), satisfies the Gromes relation [H] , 



KM + Vlir) - V^{r) = 0. 



(2) 



^'This method is often called as the Stochastic Vacuum Model. 



Notice that this relation refers both to the perturba- 
tive and nonperturbative parts of the potentials 14, (r) 
(77 = 0,1,2) and, while their perturbative parts satisfy 
this relation identically, their nonperturbative parts sat- 
isfy it in the FCM. With the definition of Vn{r) used in 
Ref. [H] the fulfilment of Eq. (O is not evident. In all 



A. M. Badalian, A. V. Nefediev, and Yu. A. Simonov 



lattice calculations (see, for example, Refs. [HHH]) re- 
lation ^ is satisfied only approximately. For a detailed 
discussion of this issue see Ref . T^ . 

The spin-dependent interactions in the form of 
Eq. ([T]) can be derived naturally in the framework of 
the FCM [17] . To this end a quark-antiquark system is 
considered and its wave function is built: 



* 



(in, out) 



{x,y\A) = *g(x)$(a;,y)*g(y), 



with nonlocal gauge invariance guaranteed by the par- 
allel transporter, 

$(a;,2;)=PexpU5 / dZf.Ap" 

Then the Green's function can be constructed as 
Gg, = {^^°;''\x,y\A)^^;f{x,y\A)),,A 

= {TvSgix, x\AMx, y)Sgiy, y\AMy, x))a, 

where Sq and Sq are the propagators of the quark and 
the antiquark, respectively, in the background gluonic 
field. Using the Feynman-Schwinger representation for 
the single-quark propagators one can see that the in- 
terquark interaction is described in terms of the Wilson 
loop W{C), with the contour C running over the quark 
trajectories, averaged over the background gluonic field. 
More specifically, the value {TiW{C)), which enters the 
quarkonium Greens' function, can be expressed through 
the correlators of the field strength tensors as 

{TtW{C)) - (riexpig f dTTf,,{z)F^,{2 



= exp^ 
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where the cluster expansion theorem [TBI (H] was used. 
The average {(■■■)) stands for connected correlators, 
for example, for the bilocal correlator, ((F(1)F(2))) = 
{F(1)F(2)) - (F(1))(F(2)}, and F^, = d^A, - d,A^ - 
ig[Afj_,Ai,] is the vacuum field strength. Obviously, 
{{F)) = {F)=0. 

The element of integration, 

dTri_,i,{z) = ds^_,v{z) - ia-j^Jdri + ia'j^JdT2, 

with aj^l = ji{jf_cju - lulfi) (i = 1, 2 for the quark and 
antiquark, respectively), contains both the element of 
the surface area and the quark spin variables accompa- 
nied by the quark proper time differentials. This is the 
most economical way to include spin-dependent inter- 
actions into consideration [181 . In the Gaussian approx- 
imation for the vacuum, when only the lowest, bilocal 



correlator is retained (see Ref. [T^ for the discussion of 
the accuracy of this approximation), one has 



{TrW) oc exp 
where 



d-K ^i,{x)dTixp{x')D ^^Xp{x - x') 
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D^,xp{x -x') = j- (TrF^,(a;)$(x, x')Fxpix'Mx' , x)) . 

^^ c 

This bilocal correlator can be expressed through only 
two gauge-invariant scalar functions, D{z) and Di{z) 
[Hi: 



Dp^\p{z) = {5p^\8i,p - S^pS^\)D{z) 



t: — (zxSiyp - Zp5\i,) 



Di{z). 



The correlator D{z) — D{zo,\z\) contains only a 
nonperturbative part and it is responsible for the QCD 
string formation at large interquark separations. The 
fundamental string tension can be expressed as a dou- 
ble integral: 



du / dXD{iy,X). 



(4) 



This correlator decreases in all directions of the Eu- 
clidean space, and this decrease is governed by the glu- 
onic correlation length T^, 

D{z) - -^ exp[-z/Tg], z = v/AM^, (5) 

where the coefficient is chosen to satisfy the relation 
(|3]) . Strictly speaking, the asymptotic form of the cor- 
relator ^ is valid at large distances, at \z\ ^Tg, while 
at \z\ < Tg the function D(z) admits a Taylor expansion 
with the parameter z^ /Tg (see, for example, Ref. |20|). 
Notice however that, for small Tg, the contribution of 
the region \z\ < Tg to the integrals defining the spin- 
dependent potentials (see Eq. ([5]) below) is negligible, 
so we stick to the simplified form ^ from the begin- 
ning. 

Another important comment concerning the corre- 
lation length Tg is the scale at which it is defined (see, 
for example, Ref. [H] for the discussion of the issue in 
relation to the Operator Product Expansion in QCD). 
By natural arguments we expect the scale of Tg to be of 
order of the average size of gluelump (average momen- 
tum), which is of the order of 1 GeV (see Ref. [22]). 

The other correlator, Z?i(z), contains both pertur- 
bative and nonperturbative contributions. Its pertur- 
bative part leads to the colour Coulomb interaction 
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between the quarks, whereas its nonpeturbative part 
was parametrised through the gluelump spectrum in 
Ref. [23]. We use this parametrisation here: 



Diiz) = 



iasCp 1 2C_Fasa'adj 



T^z 



exp[-z/T'], 



where Cf = {N^ - l)/(2iVc) - 4/3, Cadj = A^c = 3, 
(Tadj = (Cadj/Ci.)a = (9/4)a, and as = gVi^ir). Then 
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For the sake of generahty the correlation length T' is 
kept different from Tg — see Ref. [23] for the details. 

The spin-independent part of the interquark interac- 
tion follows from the area law asymptotic for the Wilson 

loop ig, 

{TrW{C)) r^expi-aSnUC)), 

where S'min is the area of the minimal surface defined 
by the contour C, 



^(C) 



dt / dpyJiww'Y-w^w'^, (6) 
Jo 



with the string profile function w^{t,(3) approximated 
by the straight-line ansatz: 



u^{t, f3) = /3xif,{t) + (1 - P)x2^{t). 



(7) 



Here xi^2(t) are the four-coordinates of the quarks at 
the ends of the string. This approximation is valid at 
least for not large excitations due to the fact that hy- 
brid excitations responsible for the string deformation 
are decoupled from mesonic excitations by the mass gap 
of order 1 GeV. 

The spin-dependent interquark interaction follows 
from the mixed terms of the form ds^j^^a\„Ff^,jF\„dT in 
Eq. dSD, 



Lso 



.(1) 



dSf^u{w)dTi(j\pD^j_^\p{w - xi) 



fl 



2), 



ds^^, = e''''daWf,{t, P)dbW^{t, l3)dtdP, a,b= {t,(3}. 

For the ansatz ([7]) one has: 

ds^i = r^dtdfi, dsik = e^kmPmdtdP, 

p^[rx {(3xi + {1 - I3)X2)], 

and thus the angular momentum enters the interaction 
through dsik- Finally, we introduce the laboratory time 
t instead of the proper quark times as dri — dt/{2^i) 
[i — 1,2), where /ii stands for the i-th quark energy 



[^ [5S]. Notice that the proper inertia of the string 
is to be taken into account when proceeding from the 
quark velocities to the angular momentum variables 
|26j . However, the effect of the string inertia and the 
deviation of the quark energy from its mass are impor- 
tant for light quarks, whereas they are negligible in case 
of heavy quarks, so that, for the heavy quarkonium, one 
has simply /i^ — rrii and Pi = rriiXi. Then, by an ex- 
plicit calculation, one can arrive at Eq. ([T]) (the details 
of the derivation can be found in Refs. [I8l[27]) with the 
following identification of the potentials Vn (r) , n = 0-4: 

Voir) = 2 / di^ dXD{X,iy)+r / dvDi{r,v), 
Jo Jo Jo 

V({r) - -2 /" diy f dxfl- - j ^(A, i^), (8) 

^2 W = - / dv I XdXD{X,iy)+r dvDi{r,v), 

T Jo Jo Jo 

Vz{r) 



d f°° 



Vi{r) 



6 / dv 

10 



D{r,v) 



1 



2 2_a_ 



Di{r,v) 



At large interquark separations, r ^ Tg,T' these yield 
simply: 



^o'(0 
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^4W = y7ras(5(3)(r). 



In particular, the static interquark potential comes out 
from Eq. ^ in the standard "linear-|-Coulomb" form. 



Vo{r)^VQQ{r) 



4 as 
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3. RESULTS AND DISCUSSION 

With the form of the potentials ([5]) we are in 
a position to fit the corresponding lattice data from 
Ref. [TT]. We have the set of four fitting parameters: 
{as,cr, Tg,r^}. In Table 1 we give the set of our fits. 
For the fits 1-3 we use the full form of the potentials 
([8]), the fit 4 demonstrates the relevance of nonpertur- 
bative interactions since only perturbative part of the 
potentials 14, (r) (n = 0-4) is retained in this case. Fi- 
nally, the data were fitted with the help of the asymp- 
totic large-distance potentials ([9]) (fit 5). Comparison 
of our fits with the lattice data is given by Figs. [T][5] (we 
use the data from Ref. [H] for the 20'^40 lattice with the 
coupling 13 = 6.0 and the lattice spacing a = 0.093 fm). 
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One can draw several conclusions from Figs. [T][5] i) 
the data clearly indicates the presence of nonperturba- 
tive contributions to the potentials even at small in- 
terquark separations, where the perturbative physics 
dominates; ii) although the simplified fit 5 approximates 
the data rather well, the fits 1-3 give a better description 
of the latter. Thus one concludes that the present data 
allow one to study the "anatomy" of the field correla- 
tors; iii) by comparing fits 1-3 with one another one can 
conclude that the data prefer small correlation lengths 
Tg and T' a good description is achieved with them 
both being 0.1 fm or below; iv) the theoretical expres- 
sions for Vn{r) (n = 0, 1, 2) given by our Eqs. ([8]) satisfy 
the Gromes relation ^ identically, while the lattice po- 
tentials, as was stressed in Ref. [H], violate it. The 
source of this violation is not clear at the moment: it 
may be a lattice artifact and disappear in the contin- 
uum limit [28j or it may be related to the contribu- 
tion of higher-order correlators taken into account dif- 
ferently in different potentials entering the Gromes re- 
lation [3]. This question deserves an additional careful 
investigation. In particular, a detailed comparison of 
the Eichten-Feinberg and Field Correlator definitions 
of these potentials is given in Ref. [IS] ; v) there is a cer- 
tain contradiction between the theoretical predictions 
and the lattice data for the potential V4(r). Our form 
of the potential V4(r) given by Eq. ([5]) is consistent with 
the delta-functional form of this potential in the limit 
of the vanishing correlation length (see Eq. ([9])). Al- 
though, for finite values of the correlation length, V4{r) 
is smeared and can become negative, the small-r be- 
haviour of our fits is essentially different from that of 
the lattice data, so this question deserves additional in- 
vestigation. 

Notice that the values of the correlations lengths ob- 
tained in the fits to the data are in good agreement with 
the predictions made for the gluelump spectrum in the 
framework of the FCM [23] . Indeed, the inverse of the 
correlation length gives the mass of the lowest gluelump 
in the corresponding channel. From the best fit 3 these 
masses constitute 3 GeV and 2 GeV for the correlators 
D{z) and Di{z), respectively. These should be con- 
fronted with the predictions 2.8 GeV and 1.7 GeV made 
in Ref. 123] (after the proper rescaling from a = 0.18 
GeV^ used in Ref. ^ to a ^ 0.22 GeV^ extracted 
from out fits). 

The fact that the vacuum correlation length is small, 
less than 0.1 fm, allows one to justify the use of the so- 
called string limit of QCD, when this length is set equal 
to zero. In this limit, the interaction of colour con- 
stituents in hadrons can be described in terms of the 
infinitely thin QCD string with the string action given 
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The set of fits for the lattice data on spin-dependent 
potentials taken from Ref. [TT]. Eqs. (|S| and (O are 
used for the fits 1-4 and fit 5, respectively. 



by Eq. ^. This approach was successfully used to de- 
scribe conventional mesons [25l |30] and baryons [31], 
hybrids 22jj glueballs [53], and gluelumps [22]. 

We conclude that the existing lattice data on the 
spin-dependent potentials in heavy quarkonia are con- 
sistent with the predictions of the FCM and that the 
gluonic correlation length extracted from these data is 
small, less that 0.1 fm. The definition and the behaviour 
of the lattice potentials Vi{r) and Vi{r) might need a 
better justification and, possibly, some improvements. 
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